Spin- and Charge Excitations of the Triangular Hubbard-Model: a FLEX-Study 
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A study of the quasi-particle excitations and spin fluctuations in the one-band Hubbard-model 
on the triangular lattice with nearest- and next-nearest-neighbor hopping is presented. Using the 
fluctuation-exchange-approximation (FLEX) results for the quasi-particle dispersion and life-time, 
the Fermi surface, and the static spin structure factor will be discussed for a wide range of dopings 
and as a function of the Coulomb correlation strength U. It is shown that the renormalization of the 
spin- and charge-dynamics is sensitive to the interplay between van Hove singularity-effects and the 
nesting, which is influenced by the next-nearest-neighbor hopping. For all dopings investigated, the 
energy-dependence of the quasi-particle life time is found to be of conventional Fermi-liquid nature. 
At intermediate correlation strength the static structure factor is strongly doping dependent, with 
a large commensurate peak at the -ftT-point for 1.35 electrons per site and weak, incommensurate 
intensities occuring at lower electron densities. The relevance of this model to the recently discovered 
cobaltates NazCo02 ■ J/H2O will be discussed. 

PACS numbers: 71.10.Fti, 71.27.+a, 74.20. Mn, 71.28. +d 
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I. INTRODUCTION 



Recent discovery of superconductivity in the cobal- 
tate Na x Co02 • J/H2O, x ~ 0.35 with a crit- 
ical temperature of T c ~ 5K— has attracted 
enormous attention to this material, stimulating a 
large number of experimenta&ii^ and theoretical 
investigationavi7i a iPi 1t )i 11 i 12 i 1i ?i 14 i 1R i 1 r. Na x Co0 2 • yH 2 is 
of space group P63/mmc, representing a hexagonal lay- 
ered structure which consists of two-dimensional (2D) 
C0O2 conduction planes, separated by H2O and Na + 
ions. The material is derived from its parent compound 
NaCo204, equivalent to a;=0.5, by simultaneous interca- 
lation of H2O and deintercalation of Na + ions. This pro- 
cess nearly doubles the c-axis lattice-constant, strongly 
enhancing the 2D character of the compound. The Co 
ions form a triangular lattice and are coordinated by 
edge-sharing oxygen-octahedra which display a shift of 
the oxygen atoms along the (1, 1, Indirection. In this 
environment the low-lying t2 fl -manifold splits into a non- 
degenerate A lg level and an E s -doublet, leaving Co 4+ in 
a low-spin, 5=1/2 state, with a single electron in the A\ g 
level while Co 3+ has a filled A% g orbitalli. In turn, x can 
be viewed as electron doping away from half-filling of a 
non-degenerate one-band model on a triangular lattice, 
with Na K Co02 • 2/H2O, x f=a 0.35 refering to an electron 
density of n « 1.35 per site. 

At present, the electronic structure of the cobaltates is 
an open issue. Yet, there is a tempting similarity to the 
high-temperature superconducting cuprates, where how- 
ever d-holes are arranged on a 2D sgware-lattice. Early 
on, the role of electron correlations in the cobaltates 
has been of interest. LDA results in an A\ g bandwidth 
of W fwl-1.4eV, while the d-shell Coulomb-correlation 
strength U is of order 5-8eV typicallj»l£. Motivated by 
this, resonating-valence-bond (RVB) 1 ^ 1 !? scenarios have 
been pursued by several researcher9§i2i±2iI4. In addi- 



tion, weak-coupling calculations, using renormalization- 
group (RG) methods have also been performed^. In this 
paper we apply an alternative method, suitable for in- 
termediate correlation strength, namely the fluctuation- 
exchange-approximation (FLEX)22i2i and investigate the 
spin- and charge-dynamics in the triangular Hubbard- 
model. In passing, we note that such analysis may also 
shed light onto the organic superconductors k-(ET) 2 X^2. 
These have been suggested to realize the Hubbard-model 
on an anisotropic triangular lattice^. 

The paper is organized as follows: in section |n] we 
discuss the model and briefly review the FLEX method. 
In section IIIII results for the renormalized dispersions, 
Fermi surfaces, the static spin structure factor, and the 
frequency dependence of the selfenergy are discussed for 
several dopings. 



II. MODEL AND METHOD 

Our subsequent analysis is focussed on the normal 
state properties of the one-band Hubbard-model on the 
triangular lattice 
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(1) 



Due to the lack of electron-hole symmetry, a choice 
has to be made to fix the sign of the the hopping ma- 
trix elements tij and the Coulomb repulsion U. Here 
we fix these parameters to model the electronic struc- 
ture of Na^OoC^ • yH 2 which we assume to be simi- 
lar to the one of the parent compound. LDA calcula- 
tions on NaCo20^ii place the Fermi energy in a CEF- 
split ti g manifold of Co d-bands with the states near the 
Fermi energy Ef to be of a\ g symmetry. Due to a weak 
unit-cell doubling a 'bilayer'-type splitting exists which 



2 




-6-4-2 CO 2 4 6 -4 



k 4 



Figure 1: a) dispersion, b), c) Fermi surface, and d) density 
of states for the thight-binding energy of eqn. for t' = 
versus t' 7^ 0. In a) [d)j horizontal[vertical] solid(dashed) lines 
refer to chemical potentials [i at n = 0.65. In b) and c) dahed 
hexagons label Brillouin zone of triangular lattice. 



equation, which is rather different from the approach pre- 
sented here. 

To study the model I|1I2I) we employ the diagrammatic 
FLEX approximation^. The FLEX is conserving in the 
sense of Kadanoff and Baym^i. It sums all particle- 
hole(particle) ladder-graphs for the generating functional 
selfconsistently and it is believed to apply up to interme- 
diate correlation strenght U/t ~ 1. The FLEX equations 
for the bare(renormalized spin- and charge-) susceptibili- 
ties X°(x s '°) I the one-particle Greens function (selfenergy) 
G(£), and the effective interaction V read 



(v m ) 



T 

N 



^ G k +q (uJn + V m ) G k (uJn) (3) 

(4) 



X\{v m ) [iTUxlivm)}- 1 
J 1 



we will neglected henceforth. LDA finds the maximum 
of the ai g -t2g bands at the T-point. Focussing at first 
on the nearest-neighbor(n.n.)-hopping t = t^, and us- 
ing an electron-language this implies that t > in eqn. 
QJ. Coulomb-correlations however, suppress double oc- 
cupancy of the a\ g -t2g level by holes. Therefore, for the 
remainder of this paper we switch to the hole-language, 
implying t > and U > 0. The corresponding tight- 
binding band and bare DOS, as well as the Fermi surface 
(FS) and chemical potential for n w 0.65 are shown in 
Fig. In the remainder of this paper all energies will be 
expressed in units of t. 

Next, it is important to realize that a dispersion with 
n.n. -hopping only is qualitatively different from the LDA. 
First, the latter exhibits a pronounced dip at the T-point. 
Second, the LDA FS for the parent compound NaCo204, 
i.e. at n=0.5, is notably hexagonal, while the tight- 
binding FS using n.n. -hopping only is practically circular 
at this doping. These facts signal a sizeable next-nearest- 
neighbor(n.n.n.) hopping t' = in the cobaltates. 

Fitting the corresponding dispersion 

e k = -2i(cos (fc x ) + 2 cos(fc x /2) cos(V3k y /2)) 

- 2t'(2 cos('Sk x /2) cos( V3fc y /2) + cos(V3fc y )) (2) 

to the LDA along the T-M direction we find t' « —0.45. 
As can be read of from Fig. ^ this does not only lead 
to the required dip in the dispersion and the hexagonal 
FS shape, but also to a signature of van Hove singular- 
ities very different from the case t' = 0. For t' = the 
bandwidth is W = 9, while for t' = —0.45 it is slightly 
increased to 9.26. The role of n.n.n. hopping has also 
been emphasized in a recent search for possible /-wave 
SC2£ using third order pertubation in U/t for the gap 



where ui n = iitT(2n + 1) and u m = ii:T2m. Eqns. l.'>l(il) 
have to be solved selfconsistently taking into account that 
G^ 1 (w n ) = [uj n — ek + /i - Sk (wn)P ■ We compute the 
Matsubara summations using the 'almost real contour' 
technique of refi^4. I.e., contour integrals are performed 
with a finite shift 17 (0 < 7 < iTir/2) into the upper 
half plane. All final results are analytical continued from 
uj + ij onto the real axis uj + i0 + by Pade approximation. 
The following results are based on FLEX solutions using 
a lattice of 64x64 sites with 4096 equidistant w-points in 
the energy range of [—30, 30]. The temperature has been 
kept at T = 0.05 . 

III. RESULTS AND DISCUSSION 

We start the discussion of our results with the quasi- 
particle dispersion Ek which is determined by — (y + 
fi — Ek(-Bk) = 0. This dispersion is shown in Fig. EI a) 
and (b) at U/W ~ 1 for t' = and —0.45 for various 
doping levels S, where S = n — 1 within the hole-picture, 
i.e. Na K Co02 • yH 2 with x ~ 0.35 corresponds to 5 m 
—0.35. As compared to Fig. ^ (a) a pronounced mass- 
enhancement of order unity at the FS crossings is clearly 
visible in Fig. [2 (a) and (b). This enhancement is due to 
low-energy spin fluctuations which are present in x q (w). 
In addition the overall shape of the dispersion is strongly 
renormalized. At large distances from the FS however, 
and due to strong quasi-particle scattering, the notion of 
the quasi-particle pole Eu may be of less relevance. 

At U = 0, t' = 0, and for 6 = +0.5 the Fermi en- 
ergy coincides with the van Hove singularity related to a 
flattening of the bare dispersion observable in Fig. (a) 
along T-M-K. Fig. [2(a) shows a shift of this flattening 
at finite U . Eg. for U = 8 we find the Fermi energy to 
coincide with the remnants of the van Hove singularity 
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Figure 2: Quasi-particle dispersion along Y-K-M-Y for 
(a) t' =0 and (b) t' — —0.45 and for various doping levels at 
(7 = 8 and T = 0.05. 

only at 6 R* +0.35. This is rather remote from the SC's 
doping level. In ref£ Fermi surface pinning of the van 
Hove singularity has also been investigated, however at 
fixed chemical potential. For t' ^ the situation is very 
different since two van Hove singularities exist in the bare 
DOS. Remnants of these remain present also at finite U. 
At t' = -0.45 and U = 8 Fig. (b) shows that for the 
SC's doping level of S « —0.35 the FS is rather close to 
the lower one of the van Hove singularities - which is ab- 
sent at t' — 0. Finally we note from Fig. (a), that for 
t' = the relative mass-enhancement at the FS along di- 
rections unrelated to the van Hove singularity, e.g. T-K 
is largest if the van Hove singularity is closest to the FS. 
For t' = -0.45 Fig. 0(b) does not show this effect. 

In Fig. the FS of the t- and t - f - U models 
are compared for S = 0, ±0.35. As compared to the 
non-interacting case at U = there are no indications 
of sizeable, correlation-induced shape-changes of the FS. 
The FS at finite and negative t' — —0.45 is hexagonally 
shaped for 8 < down to the SC's doping level. This 
is consistent with the LDA results of ref A In contrast, 
the FS at t' = evolves towards a more circular shape as 
S = —0.35 is reached. Nesting of the flat regions of the 
FS with the wave-vector Q (S) indicated in Fig. can 
enhance the static spin structure factor at Q(S). Due 
to the shape of the FS this enhancement is more pro- 
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Figure 3: Fermi surfaces for t = (thin lines) versus t' = 
—0.45 (thick lines) at U = 8 and T = 0.05 for various S. 
Outer hexagon refers to BZ. Arrow denotes a commensurate 
wave vector K = (4^, 0) with a circle labeling the half width 
of Xq at q = K at f = 0.05, (7 = 8, and f = 0. 

nounced for t' = —0.45 at 5 < 0. Interestingly, at the 
SC's doping level the nesting- vector Q is commensurate, 
i.e. Q(<J w -0.35) = K = (^,0) . 

Now, we discuss the evolution of the static spin struc- 
ture factor Re [Xk( w = 0)] with dopings S and correlation 
strength U along the path T-K-M-T in momentum space, 
contrasting the case of t' = versus t' = —0.45. Fig. 
displays the doping dependence at T = 0.05 and U — 8. 
For S = +0.35 this figure shows the static structure factor 
to be only moderately enhanced whithin a rather broad 
range of incommensurate momenta. For tf = these 
momenta are located mainly along r — M and partially 
along T — K. The zone-folded approximate nesting vector 
is located slightly off the M-Point. For t' = -0.45 and 
S = +0.35 enhanced spin fluctuations are found primarily 
along the direction T-M with Re = 0)] rather small 

within the remaining k-space. As the doping is reduced 
the maxima in Re [Xk( w = 0)] move towards the if-point 
and and develop into a sharp and commensurate peak 
at Q(S) = K as the SC's doping level of 5 = —0.35 is 
reached. The shoulder at M in Fig. El is due to the over- 
lap with the remaining peaks at symmetry equivalent po- 
sitions. Nearly all incommensurate spin fluctuations are 
suppressed. The half width of the commensurate spin 
fluctuation peak for S = —0.35 and t' = is scetched 
by the circle on the right hand side of the inner FS in 
Fig. El This circle is also a measure for the degree of 
nesting of the quasi-parallel regions of the Fermi surface. 
At 6 = —0.35, close to the SC's doping level, the com- 
mensurate peak is ^60% larger for t' = —0.45 than for 
t' = 0. At S = +0.35 however, where the nesting proper- 
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Figure 4: Doping dependence of the static spin structure 
factor Rexk(^ = 0) along the path T-K-M-T in the BZ for 
(a) t' = and (b) t' = -0.45. Notice that only at 8 = -0.35 
the commensurate peak at K is observed for large U. 



ties of the t— and t — t' — U models are rather reversed, 
the spin response is less supressed for t' = 0. 

In Fig. (a) and (b) the dependence of the static spin 
structure factor on the correlation strength U is displayed 
at the SC's approximate doping concentration for t' = 
and t' = —0.45 along a path in the BZ identical to that 
of Fig. Q] Evidently, increasing U at this doping level, 
predominantly increases Re [Xk( w = 0)] at the commen- 
surate momentum K. Additionally, one realizes that this 
increase is non-linear in U, with a rather sizeable growth 
setting in for U > 4. Finally, and similar to the dop- 
ing dependence, the overall size and sensitivity of the 
spin fluctuations to the correlation strength is larger at 
t' = —0.45 than at t' = 0. In contrast to the situation at 
S = —0.35 we have found the incommensurate spin fluc- 
tuation peaks at S — +0.35 to be enhanced only weakly 
with U increasing up to U = 8. 

Fig. shows the low-temperature frequency depen- 
dence of the imaginary part of the quasi-particle selfen- 
ergy Im£ke-Fs( w ) for 5 = 0, and ±0.35. The lattice mo- 
menta are chosen as close as possible to the FS. We find 
the selfenergy to be nearly isotropic along the FS with 
only a rather weak maximum occuring in the direction of 
the commensurate spin fluctuations. Near the Fermi en- 



Figure 5: Dependence of the static spin structure factor 
ReXktk-" = 0) on the correlation strength U along the path 
T-K-M-T in the BZ for (a) t' = and (b) t' = -0.45. 



ergy and for all dopings shown, the selfenergy is clearly 
proportional to ui 2 at low-energies which is indicative of 
normal Fermi-liquid behavior. This is in sharp contrast 
to FLEX analysis of the Hubbard-model on the square 
lattice, close to half filling. There one typically finds 
'marginal' Fermi-liquid behavior with ImEkeFS^) ^ u 
over a wide range of frequencies^*^. Therefore, along 
this line one is tempted to conclude, that the normal 
state of the cobaltates is more of conventional metallic 
nature than in the cuprates. This is even more so, if one 
realizes from Fig.|£|that the quasi-particle scattering rate 
displays its smallest curvature at the SC's doping level, 
which implies the quasi-particles to be rather well defined 
there. At larger values of S, proximity of the FS to the 
van Hove singularities enhance both, the absolute value 
of ImSkeFs(^) as well as the curvature. This effect is 
more pronounced for t' = 0. 

In summary we have investigated the one-band 
Hubbard-model on a triangular lattice with nearest- and 
next-nearest-neighbor hopping using the normal state 
version of the FLEX approximation. Adapting this 
model to the cobaltates we have shown that a sizeable 
next-nearest-neighbor hopping is necessary to describe 
their electronic structure. Over a substantial range of 
dopings which we have investigated we find a significant 
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Fermi surface mass-enhancement of order unity. This is 
due to spin fluctuations and the additional proximity of 
the Fermi surface to van Hove singularities at particu- 
lar doping levels. In contrast to the Hubbard-model on 
the square-lattice we have found the quasi-particle scat- 
tering rate to display a conventional Fermi-liquid type 
of energy-dependence. Finally we have shown that the 
static spin structure factor exhibts a large commensurate 
peak at those doping levels attributed to the supercon- 
ducting cobaltate systems. This response was found to 
be significantly enhanced by next-nearest-neighbor hop- 
ping. 
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Figure 6: Frequency dependence of — ImEk (u) near the FS 
in direction T —* K for (a) t' = and (b) t' = —0.45 and 
S = 0, ±0.35 (location of momenta k in units of K for (a): 
S = -0.35: 33/64; 6 = 0.0: 42/64; 6 = +0.35: 45/64 and (b): 
S = -0.35: 33/64; S = 0.0: 39/64; S = +0.35: 42/64 ). 
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